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Abstract 

We investigate the electric polarizabilities of neutral baryons using the method of QCD sum 
rules. The diagrams in the operator product expansion are taken into account up to dimension 6. 
We obtain different values of the polarizabilities among hyperons, the ordering of which is given 
a n > a s o > qa- The polarizability of S°, a s o, has possibly a very small value. 
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Electric and magnetic polarizabilities, labeled at and j3, are fundamental constants char- 
acterizing the hadron structure. These constants measure the ease with which a hadron 
acquires electric and magnetic dipole moments in external electromagnetic fields [Q. While 
for the proton and neutron they have long been studied through theoretically and experi- 
mentally and recent measurements have established certain values , very little is known 
about hyperon polarizabilities. However with the hyperon beams at CERN and Fermilab, 
the situation is expected to change. In particular the S polarizabilities will be soon measured 

This has induced a number of theoretical investigations. In fact predictions have been 
made in the non-relativistic quark model J7j, the heavy baryon chiral perturbation theory 
|| and the bound state soliton model ||. In the non-relativistic quark model the results are 
«s+ = 20.8, a?,- = 5.7, all in units of 10~ 4 fm 3 . The large difference in them can be traced 
back to the different quark structure of the S + and S~. While in the S + the two quark 
flavors have opposite electric charges, all valence quarks in the S~ have the same charges , 
which inhibit internal excitation. In agreement with the non-relativistic quark model, the 
heavy baryon chiral perturbation theory (HBChPT) at leading order finds «s+ > a%- ; 
namely, «s+ = 9.4(8.8) and a^- = 6.3(5.9) in the same units for the axial-vector coupling 
constants D = 0.8(0.75) and F = 0.5(0.5). Here the reason for this result is the kaon cloud 
contribution, which is strongly suppressed for the £~. In addition hyperon polarizabilities are 
smaller than the nucleon polarizabilities (a p = 13.0(12.1), a n = 11.4(10.5)) in the HBChPT. 
The large value of in the HBChPT is due to the large contribution of pion cloud. On the 
other hand the hyperon polarizabilities mainly come from the contribution of kaon cloud, 
whose mass is larger than that of pion and inhibit internal excitation under the external 
field. 

We use a completely alternative approach, the method of QCD sum rules pj| , |Tl"|j . In a 
previous work [|12|] we constructed the QCD sum rule for the electromagnetic polarizabilities 
of the neutron and found that they are related with the vacuum condensates of quark-gluon 
fields in a weak external constant electromagnetic field. We have seen that the neutron 
electric polarizability is well described with the method. In this paper we extend the method 
to hyperons and calculate the electric polarizabilities of neutral hyperons. 

We consider a two-point correlation function in a weak external constant electromagnetic 
field F^ v as follows: 



n F (p) = i J rf 4 xe^(0|T[r ?B (a;)^(0)]|0) F , (1) 

where t]b(%) is an interpolating field constructed out of quark fields such that it has the 
quantum numbers of a baryon under consideration. Following Ioffe and Smilga we 
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expand the correlation function in powers of up to order F 2 terms: 

u F ( P ) = nt°)(p) + n«( P )F^ + ng p »F^F^ + o(f 3 ). (2) 

The leading term n^°'(p) is the correlation function when the external field is absent. 11^ (jo) 
is known to be related with the magnetic moments of baryons. It is n^] (p) that is related 
with the electromagnetic polarizabilities of baryons. Requiring parity invariance, we can 
decompose US 2 } [p) into the following Lorentz structure: 



njf^a (p) = g^gvoRst {p 2 ) + zvv^n^ (p 2 " 

+g w 9uaP ! 'n-v 1 (p 2 ) + PuipgvaRv* (p 2 ) + p^Ppgvap 1 ^ (p 2 ) (3) 

+Op\P X P P gvo'H-T 1 (p 2 ) +«75WA[ff A ff n ftb 2 ) +P X Pa^P 2 (p 2 )} 

+i^ upX {g x Jn Al ( P 2 ) + P x la iiA 2 (p 2 ) + rWu,(p 2 ) + P Wn^(p 2 )]. 



Il(p 2 ) = — I dt— — + subtraction terms, (4) 



The invariant functions in the right hand side satisfy dispersion relations: 

ix Jo t — p 1 

where the subscripts S\ etc. are suppressed for brevity. One can apply the Borel transform 
B defined by 

U(p 2 ). (5) 



£[II(p 2 )] = lim v p 



9 71 • 

n— >oo 

p 2 _„ 



d(-p 2 



Here s denotes the squared Borel mass. We make a Borel transform on both sides of the 
dispersion relations: 

^ 1 r°° , 

B[U(p 2 )} = - dte- t/s lmU(t). (6) 

7T JO 

Evaluating the left hand side by an operator product expansion (OPE) in the deep Euclidean 
region we obtain the Borel sum rules. 

Let us next consider the physical content of the invariant functions in the right hand 
side of Eq.(||). The lowest intermediate state of the correlation function H(p) is expressed 
by the baryon propagator: H(p) = — \ B /\jj — Mb — S(p, F)], where Mb is the baryon 
mass, and Xb the coupling strength of the interpolating field to the baryon. S(p, F) is 
the self-energy part under the influence of the background electromagnetic field, and can 
be decomposed into possible types of Lorentz structure. The invariant functions of the self 
energy on the mass-shell (p 2 = M B ) are related to physical quantities of a baryon, such as the 
magnetic moment, the polarizabilities, via an effective Lagrangian introduced by L'vov 
This Lagrangian, which was constructed by requiring Lorentz and gauge invariance, and the 
discrete symmetries, describes the low-energy interaction of spin-1/2 composite particles with 
a soft electromagnetic field. The amplitude of low energy Compton scattering is reproduced 
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from the Lagrangian. For the details see our preceding paper Ref. [I2|. We find that the 
invariant functions and Ily are related with the polarizabilities in the following way: 



n*(f) = 


\ 2 

_A B 


n Sa (t) = 




n Vl (t) = 


X 2 


u V2 (t) = 


X 2 


u Va (t) = 


X 2 
_A B 



1 



(t - M% )3 

1 



7 [4M|Ei(0] + 



(t - M 2 ) 2 



[3M B E 1 (t)+2M|E s (t)] 



(7) 



[-8M B Ei(t)] 



(t-M 2 )3' 

' [M B E 3 (t)]l, 



(t - Ml ) 2 



[M B E 2 (t) + MiS 3 (t) + 2S y (t)] 



t-Ml 



1 



(t - M 2 )3 

1 

(* - Ml) 2 

1 



[4MlE 1 (t)] + 



(t - Ml) 2 



[E 1 (t)+2M B E s (t)] 



(t - M 2 )3 

' P 3 (t)]l. 



[4S 1 (t)-MlS 2 (t)] + 



-8Ei(t)] 



t-Ml l 



-E 2 (t) 



Ey(i). 



M, 



B 



(9) 
(10) 



(t - M 2 ) 2 



[2E 2 (t)+MlE 3 (t) + 



2Syg) 



t-M 2 



;n) 



£y and E^ on the mass-shell are related with the polarizabilities a B and /3b'- 

Ey (Ml) = (a B + /9b)/2, E 5 (M b ) = 



while Ei, £ 2 and E 3 with other quantities (T^|. Defining the following combination: 
II^V) = -2n 5l (p 2 ) -p 2 U S2 (p 2 ) + 2M B U Vl (p 2 ) + M B Tl V2 (p 2 ) +p 2 M B U V3 (p 2 ), 
we find the imaginary part of n^ as ^ to be related solely with a B : 



-1m U^ B \t) = -\%6{t - M 2 B )—. 



:i2) 



(13) 



(14) 



For the magnetic polarizabilities /3 B , it seems difficult to derive relations free from the con- 
stants which are not known well. Numerical estimation of /3 B is therefore left for a future 
work. 

We write the interpolating fields for the n, S°, E°, and A, respectively, as follows [14] : 



r)n(x) 



e abc [d al (x)C lfl d b (x)]^u c (x) 



~£abc 



s aT (x)C^s b (x) 7 57^ c (x), 



(15) 
(16) 



r) S o(x) = V2e abc {[u aT (x)C lfl s b (x)]^Yd c (x) + [d aT (x)C lfl s b (x)] l5 Yu c (x)} , (17) 
Va(x) = Jle abc {[u aT (x)C^s b (x)] l5 Yd c (x)- [d aT (x)C llM s b (x)] l5 Yu c (x)} , (18) 



where C denotes the charge conjugation and a, b and c color indices. Following Ref.JI^] the 
calculations of the Wilson coefficients are straightforward. 

Let us list the vacuum expectation values of interest and then classify them according 
to their dimensions. We are interested in the quadratic terms in F^ v . Hence, the vac- 
uum expectation value (VEV) of the lowest-dimensional operator qq is the second order 
term of the quark condensate in the external field, when we expand it in powers of F pv [J 
Four- dimensional operators are F pv F pa itself, iqSjfj,D u q and G pu G pa , where the covariant 
derivative is given by D p = d p + ig(\ a /2)A p \, and S is a symbol which makes the operators 
symmetric and traceless with respect to the Lorentz indices. The contribution of the VEV 



of the last two operators are so small that they can be safely neglected |I2]. The following 
three operators: qo pv q-F p(T , gqa pv (X a /2)G pa q, qSD p D v q have dimension 5. The contribution 
of the second operator turns out to vanish, as is for the sum rule for the masses of octet 
baryons. The Wilson coefficient of the last operator is constant and vanishes after a Borel 
transform. There are three six-dimensional operators; that is, gq r y p D u (\ a /2)G a pa , qqqq and 
QCTiiuQQ^paQ- The VEV of the first operator cannot be induced by the second-order interac- 
tion of the vacuum with an electromagnetic field, since q r y p D u (\ a /2)G a pa q has odd C-parity. 
To estimate the VEV's of the last two operators we have applied the factorization hypothe- 
sis, which is firmly based in the case of the four quark operators: (qqqq)F = 2(gg)_F • (qq)o, 
(qcr pv qqa pa q) p = {q.o , ilv q)p ■ (qa p(T q)F, where (■ ■ -)f denotes the VEV in the electromagnetic 
field and (■ ■ -)o the VEV when the external field is absent. 

For the calculation up to six- dimensional terms, we are then left with the two external- 
field-induced VEV's: {qcr py q) F) (qq) F . It is assumed in Ref. JL3| that the first of them can be 
written, in general, 

(qcr pv q) F = V4:7iaF pu xe q (qq)o, (19) 

where a is the fine structure constant and x the so-called magnetic susceptibility of quark 
condensate. We extend this hypothesis to the case of the second order in F pv and define a 
"susceptibility" as follows: 

(qq) F = (1 + AitaF pv F^<Pe 2 q )(qq) . (20) 

In order to obtain the desired sum rule, we substitute Eq.flT4|) into Eq.(H). The contribu- 
tions to the imaginary parts from higher mass states are approximated using the logarithmic 
terms in the OPE, starting at a threshold S . We calculate the left hand side of Eq.(0) by 
an OPE up to dimension 6. For hyperons, we further introduce the effect of the finite mass 
of the strange quark, m s . In the sum rule for a^o the leading terms of the mass correction 



4 In a recent work Jlq ], it is argued that (qq) f cannot be expanded in powers of F^ v for massless quarks. 
We find, however, that (qq) p can be expanded in the NJL model for massless quarks [ fl6| . In this paper we 
assume that the expansion is possible. 
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appear in the dimension 6 terms. In the sum rules for a\ and a^o we consider only the 
leading correction, which appears in the dimension 4 terms, namely m s (qq). One finally 
obtain the following results: 

4 



a 



4^3 I 



1 >iT 2 (J)el(uu} s 2 E 1 (So, s) + -vr 2 {x(6e M e d + el)}(uu) sE (S , s) 



-My 



2e u e d sE (S , s) - ^7r 4 {4</>e 2 - x 2 e 2 d }{(dd) ) 2 ] } 



--\ 2 a e~ M2js 



(21) 



a 



4vr 3l 



8n 2 (f)e 2 u (uu}oS 2 E l (So, s) + -Tc 2 {x(Qe u e s + e 2 u )}(uu) sE (S , s) 



-Ms 



32 

2e u e s sE (S , s) - y7r 4 {40e 



2 X 2 ^}((^)o) 2 ]} 



2 a " 



(22) 



{ [87r 2 0es(ss) o s 2 J Ei(5'o, s) + ^7r 2 x(3e u e s + 3e d e s + e 2 )(ss) sE (S , s 



An 3 



32 



{(e u + e d )e s - 8tt (f)m s e s (ss}o}sE (S 0} s) 

2 i „2\ 



-y7r 4 {20(e 2 + e d ) - x 2 eue d }(uu} (dd)o] } 
-A E o; s oe , 



(23) 



a 



12tt 3 



{[87T 2 0{2e 2 (MM) o + 2e 2 (Jrf) - e 2 (ss) }s 2 £i(.So, s) 



-M A 



4 

-7r 2 x{2(e n + 3e d + 3e s )e u (uu) + 2(e d + 3e u + 3e s )e d (dd) 
-(e s + 3e u + 3e d )e s (ss) }s J E (S'o, s) 

{e s e u + e s e d + 4e u e d + M 2 (j)m s {2e 2 u (uu) + 2e 2 d (dd) - 3e 2 s (ss) )}sE (S , s) 
64 
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32 



vr 4 0{2(e 2 + e 2 s )(uu) (ss) + 2(ej + e 2 s )(dd} (ss} - (ej + e 2 u )(dd} (uu} } 



+—7i 4 x 2 {2e u e s (uu) (ss) + 2e d e s (dd) (ss) - e u e d (uu) Q (dd) }} } 



--X 2 A a A e 



■Ml/s 



(24) 



where E n (S ,s) = 1 - exp{-S /s) Ek=o( s o/ s ) k / kl lt is wel1 known that from U^(p) in 
Eq.(Q) one can get the Borel sum rules for \ 2 B listed below. 

s 3 E 2 (S , s) + ^(^G^osEoiSo, s) + %\(dd) ) 2 = 2(2 7 r) 4 A 2 v e- M -/ s , (25) 

7T O 

s 3 E 2 (S , s) + 7r 2 (^G 2 ) sE (S , s) + ^\{dd) Q ) 2 = 2(2^X1^'% (26) 

7T 6 
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s 3 E 2 (S , s) + 8n 2 m s (ss) sE (S , s )+n 2 (-G 2 ) sE (S , s ) + %\(uu) ) 2 = 2(2vr) 4 A|e- M ^, 

no 

(27) 

s 3 E 2 (S , s) - lrr 2 m s [2((uu) + (dd) ) - 3(ss) ]sE (S , s) + n 2 (^G 2 ) sE (S , s) 

3 7T 

+ ^vr 4 [2((^) + (dd) )(ss) - <t*«>o<&i>o] = 2(2ny\ 2 A e- M '/ s , (28) 

where ((a s /7r)G 2 ) denotes the gluon condensate. 

Here, let us compare the sum rules for og, Eqs.(^T])~(^p. In the n and 5° the left-hand 
sides are similar to one another because of the same isospin-1/2 structure of the interpolating 
fields in Eqs.([T5|) and ([0J), while the interpolating fields of Eqs. QT7D and (13) have isospin- 
1 and 0, respectively and the above similarity is not seen. In the limit of flavor SU(3) 
symmetry, the n and S° exactly coincide with each other, because of e s = e^. Furthermore, 
up to dimension 5 we find the generic results: 

4_ 

Ot-n = a E° — o«A = 4«so. (29) 
o 

It is interesting to know that not only the leading terms but also higher dimensional terms 
satisfy the above relation. 

Let us now estimate olb numerically. We eliminate in Eqs.(pTD^(p4|) by using Eqs.(p5|)^(p8| 
For the susceptibilities we use the values evaluated in Ref . [|12jl : 

X = -2.58 GeV~ 2 , = 1.66 GeV" 4 . (30) 

We note, however, that there is an uncertainty due to the variation of the susceptibilities 
with the change of s: In the reliable interval of s, x varies from -2.76 GeV~ 2 to -2.44 GeV~ 2 
,6 from 1.49 GeV -4 to 1.91 GeV -4 . The vacuum condensates of the operators and the mass 



of the strange quark are taken from Ref. [0 , and are 

{uu) = (dd) = (-225 MeV) 3 , (ss) = 0.8{uu) , {—G 2 ) = (360 MeV) 4 , m s = 150 MeV. 

7T 

In Fig.l we show the squared Borel mass s dependence of olb- We see that, for a n , a=o and 
a a, the plateau develops for the suitable values of the continuum threshold parameter So. 
This implies that the sum rules for a n , a s o and a a work very well and are reliable. Indeed, 
we see that dominant contributions in the OPE come from the leading terms and the higher 
dimensional terms are suppressed. The convergence of the OPE hence appears to be good. 
Taking the maximum value of the curves, we predict that 

a n = 14.0 x 10" 4 fm 3 , (31) 
ct H o = 10.8 x 10~ 4 fm 3 , (32) 
a A = 6.2 x 10~ 4 fm 3 . (33) 
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These results have uncertainties arising from the errors in the susceptibilities. With the 
variation of the susceptibilities, a n varies from 12.1 to 16.8, a^o 9.1 to 13.5, and 6t\ 5.0 to 
8.0. 

On the other hand, the Borel curve of a^o are not stabilized for any choice of So due 
to apparently slow convergence of the OPE. Indeed, the dimension 6 term is significantly 
large and contributes to a^o negatively. In view of this, no definite conclusion for a-^o can 
be drawn. However, we can say that a^o might take a small value. 

It is interesting to compare our results with those obtained in the heavy baryon chiral 
perturbation theory (HBChPT) ||. The QCD sum rules give smaller values for the electric 
polarizabilities of hyperons as compared with that of neutron. This is in agreement with 
the HBChPT prediction. In the HBChPT, as already mentioned, the splitting is due to 
the contribution of kaon cloud. It can be seen that the ordering of magnitudes for strange 
baryons is reversal in the two calculations; namely, in the HBChPT we find a^o > a\ > a^,o, 
while a=o > > a^o in the QCD sum rules. It is known that the effect of the decouplet 
states is large, so that it is highly desired to obtain the result of the HBChPT with including 
the decouplet states fl~7H . 

For charged baryons, there newly appear in the phenomenological side a triple pole term 
which is proportional to squared charge of the baryon and involves large factor. Within the 
calculation up to dimension 6 in the OPE, this term makes the Borel stability worse. We 
therefore need to evaluate the higher dimensional terms. This task is presently under way. 

In summary, we have presented a calculation of the electric polarizabilities of neutral 
baryons in the framework of the QCD sum rule. The operators up to dimension 6 were 
retained in the operator product expansion. The sum rule predicts that the ordering of the 
polarizabilities is a n > a=o > 6t\ > o^o, and implies that a^p has a possibility to be very 
small. In particular, we find that in the flavor SU(3) symmetric limit, a n = a=o = 4o:a/3 = 
4aso, up to dimension 5 terms in the OPE. The calculations performed in the HBChPT 
have led still different predictions. However, the future experimental data from CERN and 
Fermilab could be of help to discriminate the different approaches. 

One of the authors (Y.K.) would like to thank O. Morimatsu, M. Oka and D. Jido for 
helpful discussions. Y.K. thanks also the members of the theory group of Institute of Particle 
and Nuclear Studies for their hospitality during his stay there. 
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Figure Caption 

Fig.l 

The squared Borel mass s dependence of olb- The values of the continuum threshold param- 
eter Sq are taken to be 3.6 GeV 2 for the n, 4.2 GeV 2 for the S°, 6.7 GeV 2 for the A and 
oo GeV 2 for the S°. 
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